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In  the  design  of  the  direct  methanol  fuel  cell  and  the  evaluation  of  new  materials  and  their  appropri¬ 
ateness  for  inclusion,  it  is  helpful  to  consider  the  impact  of  material  properties  on  the  performance  of  a 
complete  system:  to  some  degree,  methanol  crossover  losses  can  be  mitigated  by  proper  system  design. 
As  such,  an  analytical  model  is  developed  to  evaluate  the  methanol  concentration  profile  across  the 
anode  backing  layer  and  membrane  of  the  direct  methanol  fuel  cell.  The  model  is  integrated  down  the 
anode  flow  channel  to  determine  fuel  utilization  as  a  function  of  the  feed  concentration,  backing  layer 
properties,  and  membrane  properties.  A  minimum  stoichiometric  ratio  is  determined  based  on  main¬ 
taining  zero-order  methanol  kinetics,  which  allows  the  fuel  efficiency  to  be  optimized  by  controlling 
these  physical  properties.  This  analysis  is  then  used  to  estimate  the  required  flow  rates  and  the  size  of 
system  components  such  as  the  methanol  storage  tank,  based  on  the  minimum  methanol  flow  rate  that 
those  components  must  produce  to  deliver  a  specified  current;  in  this  way,  the  system-level  benefits  of 
reduced  membrane  crossover  can  be  evaluated. 

©  201 1  Published  by  Elsevier  B.V. 


1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  an  electrochem¬ 
ical  energy  conversion  device  that  operates  by  converting 
methanol  and  water  to  carbon  dioxide  and  protons  at  the  anode, 
CH3OH  +  H20 ±5  C02  +  6H+  +  6e_,  and  by  reducing  oxygen  at  the 
cathode,  02  +4H+  +  4e-  ^2H20.  The  primary  advantage  of  the 
DMFC  is  that  it  can  be  operated  with  a  liquid  feed,  and  there¬ 
fore  does  not  require  a  reformer  or  a  bulky  gas  storage  tank.  Also, 
because  the  oxidation  reaction  produces  six  electrons  per  methanol 
molecule,  the  DMFC  delivers  a  high  specific  energy.  The  simplic¬ 
ity  and  compactness  of  the  fuel  delivery  system  makes  the  DMFC 
attractive  for  portable  applications. 

A  paramount  problem  in  the  implementation  of  the  DMFC  is 
methanol  crossing  over  to  the  cathode,  where  it  is  oxidized,  result¬ 
ing  in  a  depolarization  of  the  positive  electrode  and  a  waste  of  fuel. 
The  origin  of  crossover  lies  in  the  fact  that  methanol  is  completely 
soluble  in  water  and  that  water  is  needed  to  swell  the  proton- 
exchange  membrane  and  to  impart  protonic  conductivity  across 
the  membrane.  One  way  to  mitigate  this  problem  is  to  create  a  low- 
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crossover  membrane  to  replace  standard  perflourosulfonic  acid 
membranes  such  as  Nation®.  Increasing  the  efficiency  of  the  anode 
reaction  through  the  design  of  the  catalyst  layer  can  also  minimize 
crossover.  Finally,  optimizing  the  design  of  the  anode  flow  chan¬ 
nel  and  operating  parameters  such  as  methanol  feed  concentration 
and  temperature  is  crucial. 

During  the  past  two  decades,  numerous  approaches  to  model¬ 
ing  the  DMFC  have  been  taken.  Most  of  these  models  are  numerical 
or  semi-empirical,  with  only  a  handful  being  analytic,  and  most 
focus  on  polarization,  power  output  and/or  kinetics.  Among  the 
analytical  models,  the  majority  are  one-dimensional,  single-phase 
models  [1-14]  with  a  few  multi-dimensional  models  [15,16].  There 
are  also  many  DMFC  models  that  consider  explicitly  the  engineer¬ 
ing  problems  of  fuel  utilization  and  efficiency  [17-27],  but  none  of 
these  are  analytic  in  nature.  There  are  also  numerous  models  in  the 
literature  that  consider  crossover,  including  a  majority  of  the  ana¬ 
lytic  models  [1-5,9,10,12,13,15],  but  these  models  do  not  consider 
the  system-level  problem  of  fuel  utilization  and  its  relationship  to 
crossover. 

In  this  paper,  we  seek  to  calculate  the  magnitude  of  methanol 
crossover  for  different  membrane  and  backing  layer  properties 
when  we  design  the  system  around  optimized  fuel  efficiency.  First, 
we  develop  analytical  solutions  for  the  methanol  concentration 
in  the  anode  flow  channel,  anode  backing  layer,  and  membrane 
based  on  the  methanol  feed  concentration  and  the  methanol  mass 
transfer  coefficients  in  the  backing  layer  and  membrane.  We  then 
set  a  minimum  allowable  stoichiometric  ratio  in  order  to  ensure 
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diffusion  diffusion 

only  and  drag  no  methanol 


Fig.  1.  Cross-section  of  the  cell  including  a  qualitative  rendering  of  the  methanol 
concentration  profile. 


Fig.  2.  Membrane  concentration  profile  for  different  values  of  phi.  The  membrane 
concentration  is  normalized  according  to  the  anode  concentration,  and  the  mem¬ 
brane  position  is  also  normalized. 


that  the  system  is  always  operated  with  sufficiently  high  methanol 
concentrations  so  that  the  anode  reaction  remains  zero-order  with 
methanol  concentration.  The  inverse  of  the  stoichiometric  ratio  is 
the  fuel  efficiency,  so  by  setting  a  minimum  value  on  the  amount 
of  fuel  delivered,  we  set  a  maximum  value  on  the  fuel  efficiency. 
We  are  then  able  to  optimize  the  membrane  and  backing  layer 
properties  to  deliver  the  maximum  possible  fuel  efficiency.  Finally, 
we  demonstrate  how  to  determine  the  size  of  system  components, 
especially  the  methanol  storage  tank,  for  a  cell  operating  at  its  mini¬ 
mum  stoichiometric  ratio.  The  sizes  of  these  components  scale  with 
the  total  flow  rate  of  methanol  into  the  cell,  so  by  calculating  the 
minimum  flow  rate  corresponding  to  the  maximum  efficiency,  we 
can  estimate  their  minimum  sizes. 

It  is  important  to  note  that  this  model  is  concerned  only  with 
methanol  concentration  and  does  not  employ  electrochemical 
equations  to  calculate  polarization  or  power  density.  Therefore,  we 
do  not  consider  proton  transport  or  the  effects  of  ohmic  losses  in  the 
membrane.  Any  effects  stemming  from  water  transport  or  carbon 
dioxide  gas  removal  on  the  anode  side  are  ignored  or  incorporated 
into  an  overall  effective  diffusion  coefficient.  The  cathode  is  not  con¬ 
sidered  explicitly  because  any  methanol  that  reaches  the  cathode 
catalyst  layer  is  assumed  to  be  immediately  oxidized.  The  methanol 
is  assumed  to  be  an  ideal  incompressible  liquid  in  an  isothermal 
environment  such  that  any  pressure  or  temperature  effects  can  also 
be  ignored.  Much  work  has  been  done,  most  notably  by  Scott  et  al. 
[28]  and  Wang  and  Wang  [29],  to  show  that  two-phase  effects  can 
be  important  even  in  primarily  liquid-feed  DMFCs,  but  this  model, 
like  previous  analytical  models,  does  not  consider  those  effects. 
Instead,  we  assume  that  the  behavior  of  the  backing  layer  can  be 
characterized  simply  by  an  effective  diffusion  coefficient  and  by  its 
thickness. 

2.  Methanol  concentration 

2.1.  Anode  backing  layer 

The  methanol  concentration  in  the  cell  can  be  broken  down  into 
two  orthogonal  concentration  profiles:  the  concentration  profile 
across  the  thickness  of  the  cell  (Figs.  1  and  2)  and  the  concen¬ 
tration  profile  down  the  flow  channel  as  methanol  is  consumed 
through  reaction  and  crossover  (Figs.  3  and  4).  To  find  the  con¬ 
centration  profile  in  the  anode  backing  layer,  we  first  assume  that 
any  methanol  entering  the  anode  backing  layer  is  either  consumed 
in  the  desired  reaction  at  the  anode,  thereby  producing  six  elec¬ 


trons,  or  crosses  over  into  the  membrane.  In  the  development  of  the 
model  that  follows,  we  assume  that  all  fluxes  and  current  densities 
are  defined  per  unit  of  superficial  area,  and  transport  properties 
are  effective  transport  properties,  modified  from  bulk  values  by 
appropriate  corrections  for  volume  fraction  and  tortuosity.  We  also 
neglect  convection  across  the  thickness  of  the  cell.  Assuming  the 

anode  cathode 

backing  backing 


crosses  over  to  the  cathode. 


Fig.  3.  Methanol  consumption  in  the  fuel  cell. 


Fig.  4.  Normalized  feed  concentration  as  a  function  of  normalized  distance  down 
the  channel  with  v  =  2  and  \fs  =  5. 
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only  transport  mechanism  at  work  in  the  backing  layer  is  diffusion, 
the  methanol  flux  density  is 

^backing  _  *  ,  h  _  D backing  ^  ^ 

Me0H  -  6F  +  6F  “  Lbacking  {Cfeed  Canode) '  11  j 

The  flux  is  constant  at  steady  state,  so  V  •  N^aeCQ  =  0,  which  implies 
that  the  resulting  concentration  profile  in  the  backing  layer  is  linear. 

When  the  concentration  at  the  anode  goes  to  zero,  the  crossover 
current  disappears,  and  the  limiting  current  in  the  backing  layer  is 

ihm  =  6F(^l)cfeed.  (2) 

y  Lbacking  J 

Therefore,  the  maximum  rate  of  transport  from  the  flow  channel  to 
the  electrode  and  membrane  is  set  by  Cfeed ,  which  is  the  methanol 
concentration  in  the  flow  channel,  and  the  backing  layer  properties 
D backing  and  Lbacking.  The  crossover  current  can  now  be  described  as 
a  function  of  this  limiting  current 

i  _  czr  I  D backing  \  (  ^  ^  \  4  _i  I  Cfeed  ~  C anode  \ 

lx  —  Or  \  I  (C/eed  ^anode )  *  —  *lim  I  )  *• 

\Lbacking  J  \  cfeed  J 

(3) 


2.2.  Membrane 


We  assume  that  the  concentration  of  methanol  at  the  cathode 
is  driven  to  zero,  and  define  the  current  as  positive  (from  the  anode 
toward  the  cathode)  because  the  oxidation  reaction  is  occurring 
at  the  anode.  In  order  to  keep  this  model  as  simple  as  possible, 
we  choose  to  ignore  variations  across  the  thickness  of  the  anode 
catalyst  layer  when  considering  the  methanol  concentration.  We 
consider  this  assumption  valid  in  this  model  because  the  catalyst 
layer  is  approximately  an  order  of  magnitude  thinner  than  both  the 
membrane  and  backing  layer.  Previous  work  [4-6,9]  justifies  this 
approach,  as  these  papers  show  that,  even  at  high  current  densities, 
the  change  in  methanol  concentration  across  a  thin  catalyst  layer  is 
much  smaller  than  it  is  across  either  the  membrane  or  the  backing 
layer. 

The  methanol  flux  across  the  membrane  is  driven  by  diffusion 
and  electroosmotic  drag.  Both  forces  will  drive  methanol  toward 
the  cathode,  and  the  methanol  flux  density  is 

iVMeOH  -  ~Umemdx  +  p  *  W 

We  assume  that  the  electroosmotic  drag  of  methanol  is  directly 
proportional  to  the  methanol  concentration,  §Me0H  =  ?MeoH  ' c; 
with  these  assumptions  incorporated  into  the  model,  we  can  solve 
for  the  concentration  profile  in  the  membrane.  The  constant  ^e0H 
is  equivalent  to  the  factor  £h2o/ch2o>  where  §h2o  is  the  drag  coef¬ 
ficient  for  pure  water  and  cH2o  is  the  water  concentration  in  a 
membrane  exposed  to  pure  liquid  water.  We  assign  values  of 
§h2o  =  2.5  [30]  and  Ch2o  =  0.0546  mol  cm-3  [31  ]  for  pure  water 
at  60  °C  -  neglecting  the  effect  of  changing  methanol  concentration 
-  so  that  ¥Me0H  =  46  cm3  mol-1.  The  flux  of  methanol  across  the 
thickness  of  the  membrane  must  be  uniform  at  steady  state,  which 
implies  that  V  •  =  0.  Therefore,  the  concentration  profile  is 

defined  by 

d2C  ^MeOH  dc  _  q 
dx2  DmemF  dx 

which  after  applying  the  boundary  conditions  c  =  canode  at  x  =  0  and 
c  =  0  at  x  =  Lmem,  has  the  solution 


Canode 

1  -exp[-0] 


1 1  -  exp 


x 

Lmem 


(6) 


0  =  ( Lmem/DmemXtij'Me0H/F )  is  a  dimensionless  ratio  of  the  flux  due 
to  exclusively  to  electroosmotic  drag  at  a  reference  concentration 


relative  to  the  limiting  value  of  the  flux  due  to  diffusion  alone.  As 
the  plot  in  Fig.  2  shows,  increasing  the  strength  of  the  drag  force 
relative  to  the  diffusion  force  increases  the  amount  of  methanol  in 
the  membrane  and  stretches  the  concentration  profile  toward  the 
cathode. 

We  have  information  about  the  rate  at  which  methanol  can 
move  across  the  backing  layer,  and  the  rate  at  which  any  methanol 
that  crosses  the  anode  catalyst  layer  will  move  across  the  mem¬ 
brane,  but  we  seek  to  describe  these  rates  as  a  function  of  the  two 
layers’  transport  and  geometric  properties,  as  well  as  the  local  value 
of  the  concentration  in  the  flow  channel.  To  simplify  the  mathemat¬ 
ics  from  here,  we  define  another  dimensionless  parameter  n,  which 
is  a  function  of  and  the  constants  Lmem  and  Dmem : 

n  =  LmemN =  0  exp  (0) 

F^memC  anode  exp  (0)  —  1 

n  is  a  dimensionless  ratio  of  the  flux  across  the  membrane  subject 
to  the  operating  conditions  of  the  system  relative  to  the  value  of 
the  flux  that  would  be  present  without  electroosmotic  drag.  Now 
Canode  can  be  defined  in  terms  of  the  cross  over  current  ix: 


Substituting  canode  back  into  Eq.  (3),  we  can  solve  for  ix  in  a  different 
form 


=  (9) 

where  u  =  (Dmem/Lmem)(0exp(0)/exp(</>)- 'l){LbackingIDbacking)  is  a 

dimensionless  ratio  of  the  limiting  current  in  the  membrane  rela¬ 
tive  to  the  limiting  current  in  the  backing  layer,  v  serves  as  a  proxy 
for  the  amount  of  methanol  that  crosses  over  to  the  cathode. 

For  completeness,  it  should  also  be  noted  that  substituting  Eq. 
(4)  for  ix  into  Eq.  (9)  and  solving  for  canode  produces  a  simple  rela¬ 
tionship  between  canode  and  Cfeed. 


Can“te-(1  ium)  (l+y)9^' 


(10) 


Further  substitution  of  Eq.  (2)  for  iUm  results  the  following  expres¬ 
sion  for  Canode  • 


_  ^ cfeed  c/eed 

Canode  -  W+v) 


(11) 


This  relationship  is  not  used  further  in  this  paper,  but  it  does  show 
that  Canode  dependent  on  Cfeed  at  the  corresponding  location  in 
the  flow  channel,  the  initial  feed  concentration  cj?ed,  the  limiting 
current  at  the  anode  entrance  normalized  by  the  desired  current 
0  =  and  the  dimensionless  ratio  v. 


2.3.  Anode  flow  channel 


Assuming  the  velocity  v  of  methanol  in  the  flow  channel  is  con¬ 
stant,  the  total  flow  rate  of  methanol  in  the  anode  flow  channel  is 
given  by 

=  cfeed  Channel  V.  (12) 

and  the  feed  concentration  down  the  channel  is  given  by 

dCfeed  — 

_ ^backing  f  ^J-channef  \  _  _  f  +  be  \  f  ^J^channel\ 

dz,  ~  MeOH  \AchannelVJ~  U  F  J  \Achannel  V)  ’ 

(13) 

where  2!  =  zlLchannei  is  the  distance  down  the  channel,  and  w  is  the 
channel  separation  width,  that  is,  the  width  of  electrode  perpen¬ 
dicular  to  the  channel  flow  that  is  fed  by  a  single  channel.  wLchannei 
is  the  area  of  the  anode  that  receives  the  methanol  flowing  through 
the  channel.  The  effects  of  other  parallel  channels  are  ignored. 
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While  we  allow  for  the  crossover  current  ix  to  vary  with  position, 
we  make  the  simplifying  assumption  that  current  density  i  is  inde¬ 
pendent  of  position  along  the  anode  flow  field  channel.  While  not 
strictly  accurate,  a  cross-flow  anode  and  cathode  configuration  and 
zero-order  methanol  kinetics  suggest  that  local  current  densities 
are  not  a  function  of  local  methanol  concentration  and  therefore 
might  be  treated  as  independent  of  position  in  the  anode  chan¬ 
nel.  This  significantly  simplifies  the  analysis  and  offers  simplicity 
in  examining  the  overall  performance  of  the  system. 

After  substituting  Eq.  (9)  for  ix  and  applying  the  boundary  con¬ 
dition  cfeed°  at  z  =  0  with  =  zfjm  /  i,  the  solution  to  Eq.  (12)  is 

9eed  =  [c%ei  +  ^  -  Cff-  (14) 


Table  1 

Physical  properties. 


Parameter  (unit) 

Value 

Reference 

9  (dimensionless) 

10“4 

Assumed 

£>MeOH-H20  (cm2  S_1 ) 

IQ-1.416  -999.778/T 

Yaws  [34] 

e  (dimensionless) 

0.7 

Wang  and  Wang  [29] 

^backing  (cm2  S  ^  ) 

t^MeOH-H2o  £  15  (base 
value) 

Wang  and  Wang  [29] 

Dmem  (cm2  S”1 ) 

4.9  x  10-6  (base 
value) 

Kauranen  and  Skou  [33] 

t backing  (cm) 

0.0180  (base  value) 

Meyers  and  Newman  [32] 

tmem  (cm) 

0.0179 

Meyers  and  Newman  [32] 

Ch2o  (mol cm-3) 

0.0546 

CRC  Handbook  [31] 

£h2o  (dimensionless) 

2.5 

Zawodzinski  et  al.  [30] 

^MeOH  (dimensionless) 

£h2o/Ch20 

Calculated 

A  more  useful  formulation  of  the  feed  concentration  is  to  define 


As  a  check  on  the  reasonableness  of  this  solution,  we  can  find  co  in 
the  limit  where  the  crossover  (in  this  case  represented  by  u)  goes 
to  zero: 


limco  =  1  -  — .  (16) 

S 


As  expected,  co  is  unity  at  z!  =  0  and  goes  to  zero  at  the  end  of  the 
channel.  The  parameter  s  in  the  equations  above  is  defined  as 


[  Seed  \ 

f  A channel  V  \ 

{ i/6F  1 

V  L channel  w  / 

(17) 


s  is  the  stoichiometric  ratio  for  the  methanol  flow:  the  rate  at  which 
methanol  is  delivered  to  the  inlet  of  the  cell  relative  to  the  amount 
of  methanol  that  is  consumed  in  the  preferred  electrochemical 
reaction.  Hence,  it  represents  the  inverse  of  the  fuel  efficiency. 
This  ratio  is  the  critical  point  in  this  analysis,  and  in  the  following 
sections  we  use  it  to  optimize  the  system-level  fuel  efficiency. 


where  6  =  Ciim/c°ed,  which  we  set  equal  to  10-4  throughout  this 
study.  To  find  the  maximum  fuel  efficiency  that  will  accommodate 
this  minimum  feed  concentration,  we  substitute  Eq.  (17)  into  Eq. 
(14),  set  z!  =  1  so  that  we  are  at  the  end  of  the  channel,  and  solve  for 
s.  The  result  is 


(1  +  u)ln  [(1  +  lm/0/(1  +  v)  (l  +  0Ln/f)] 

Placing  this  lower  limit  on  s  places  an  upper  limit  on  fuel  efficiency 
because  efficiency  (utilization)  is  1/s. 

In  order  for  s  to  be  physically  meaningful,  it  must  be  positive 
because  all  of  the  parameters  that  comprise  it  in  Eq.  (16)  are  pos¬ 
itive  definite.  Therefore,  the  term  inside  the  natural  logarithm  in 
the  denominator  must  be  greater  than  unity,  which  requires  that 
\[f  >  1  /( 1  -  0  -  v  ^  0).  This  result  makes  sense  if  we  reverse  it  to  read 
0  +  u0  +  l/^<l.  The  term  on  the  left  is  equal  to  co  as  in  Eq.  (17). 
Because  co  is  the  ratio  of  the  feed  concentration  in  the  channel  to 
the  initial  feed  concentration,  it  must  be  less  than  unity.  What  we 
find  is  that  for  low  feed  concentrations,  the  current  must  be  kept 
low  or  the  backing  layer  must  be  very  thin  and  porous.  Otherwise, 
the  anode  concentration  will  not  be  kept  above  its  minimum  value 
(Table  1). 


3.  Fuel  efficiency 


3.2.  Optimization  of  backing  layer  properties  for  varying  current 
loads 


3.1.  Upper  limit  on  fuel  efficiency 


In  a  region  with  sufficiently  low  methanol  concentration  at  the 
anode,  the  anode  potential  must  increase  sharply  in  order  to  sus¬ 
tain  a  current  as  the  methanol  oxidation  reaction  transitions  from 
zero-order  to  first-order  [31].  Rather  than  model  the  complicated 
transition  to  first-order  methanol  kinetics,  which  would  require  a 
substantial  overpotential  at  the  anode  and  a  corresponding  loss  in 
performance,  we  ensure  that  the  optimized  system  remains  above 
any  such  transition  by  requiring  the  feed  concentration  (repre¬ 
sented  by  the  dimensionless  concentration  co)  to  remain  above 
a  minimum  value  throughout  the  flow  channel.  In  doing  so,  we 
impose  a  constraint  on  the  minimum  allowable  stoichiometric  ratio 
s,  which  in  turn  constrains  the  values  of  v  and  xj/.  Fig.  4  reveals  that  if 
the  stoichiometric  ratio  is  too  small,  there  will  be  no  methanol  in  the 
flow  channel  at  the  stack  exit  (z!  =  1 ).  This  situation  occurs  when  the 
system  demands  uniform  current  density  beyond  the  point  where 
the  flow  channel  can  supply  methanol  to  the  anode  at  a  sufficient 
rate. 

Returning  to  Eq.  (3)  relating  i  and  ix  to  Cfeed  and  canode ,  we  sub¬ 
stitute  Eq.  (9)  for  ix  and  cUm  for  canode  to  find  the  feed  concentration 
in  terms  of  the  minimum  concentration  c/im : 


Cfeed 

C0 

Lfeed 


^-(1  +  v)  +  \  or®  =  0(1  +  v)  + 

feed  * 


1 


(18) 


Now  that  we  have  a  method  for  finding  the  maximum  fuel 
efficiency  under  a  given  set  of  conditions,  we  can  look  to  opti¬ 
mize  the  efficiency  for  the  variety  of  current  loads  that  the  fuel 
cell  must  support.  If  fuel  efficiency  is  the  only  concern,  then  the 
best  choice  is  to  lower  the  membrane  permeability  as  much  as 
possible  and  always  pass  a  small  current.  However,  the  cell  is 
always  subject  to  a  load  requirement,  and  often  that  requirement 
varies. 

In  the  following  demonstration,  we  first  optimize  the  fuel  effi¬ 
ciency  for  a  standard  0.05  A  cm-2  load  by  varying  the  limiting 
current  of  the  backing  layer.  Then  we  adjust  the  limiting  cur¬ 
rent  based  on  the  design  constraint  that  the  cell  must  be  able  to 
deliver  high  current  (0.5  A  cm-2  here)  with  reasonable  efficiency 
while  still  staying  close  to  the  optimum  efficiency  at  the  standard 
load. 

The  width  of  the  membrane  is  taken  to  be  0.0179  cm  as  in  dry 
Nation®  1 1 7  [32  ].  As  a  reference  value  for  the  diffusion  coefficient  in 
the  membrane,  and  we  refer  to  the  superficial  diffusivity  measured 
by  Kauranen  and  Skou  [33]: 

Dmem  =4.9  x  10-6cm2  s_1  at60°C  (20) 

Our  reference  point  for  the  effective  diffusion  coefficient  in  the 
backing  layer  is  taken  from  Wang  and  Wang  [29],  using  their  mea¬ 
sured  porosity  of  £  =  0.7  and  the  diffusion  coefficient  of  methanol  in 
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water  DMe0H-H2o  =  1 0_1  ■416-"9-778/T  Cm2  s”1  measured  by  Yaws 
[34]: 

D backing  ~  ^Me0H-H20  £  15  =  10-1 -416-999.778/7^ .7)15cm2  S-1 

or  2.242  x  10“5cm2  s”1  at60°C.  (21) 

In  each  case,  we  fix  the  feed  concentration  at  the  channel  input  and 
vary  the  limiting  current  from  just  below  the  desired  load  value 
(0.05  A  cm-2  and  0.5  A  cm-2)  to  twenty  times  that  value.  We  then 
plot  the  maximum  fuel  efficiency  as  a  function  of  the  limiting  cur¬ 
rent  for  a  variety  of  feed  concentrations.  Referring  to  Eq.  (2 ),  the  lim¬ 
iting  current  is  directly  proportional  to  the  mass  transfer  coefficient 
D  backing  I  hacking  when  the  initial  feed  concentration  is  fixed.  Previ¬ 
ous  experimental  work  [3]  has  shown  that  the  mass  transfer  coef¬ 
ficient  can  change  with  the  methanol  feed  concentration  in  a  two- 
phase  regime,  but  the  effect  is  small  enough  that  we  neglect  it  here. 

Optimizing  the  fuel  efficiency  requires  delivering  enough 
methanol  to  the  anode  to  produce  the  desired  current  without 
delivering  so  much  that  crossover  to  the  cathode  becomes  a  sig¬ 
nificant  factor.  Hence,  the  backing  layer  properties,  which  define 
the  limiting  current  for  a  fixed  methanol  concentration,  must 
be  tuned  to  the  membrane  properties,  the  desired  load,  and  the 
input  methanol  concentration.  In  both  the  high-current  and  the 
low-current  cases,  we  can  allow  the  limiting  current  to  approach 
the  load  value,  but  the  model  will  give  nonphysical  results  if  we 
demand  a  load  higher  than  the  limiting  current. 

The  efficiency  increases  sharply  until  the  limiting  current 
reaches  about  four  times  the  load  value  (0.2  A  cm-2  in  the  low  cur¬ 
rent  case  and  2  A  cm-2  in  the  high  current  case)  and  then  levels  off 
or  in  some  cases  declines.  In  this  scenario,  operating  the  cell  at  a 
low  feed  concentration  such  as  0.1  M  or  0.5  M  with  a  limiting  cur¬ 
rent  at  the  anode  of  about  2  A  cm-2  is  optimal.  The  cell  would  run 
at  its  maximum  efficiency  during  normal,  low  current  operation 
while  still  maintaining  reasonable  efficiency  when  a  higher  load  is 
demanded  (Fig.  5) 

3.3.  Optimization  of  membrane  properties 

Now  we  can  apply  the  same  method  to  study  how  fuel  efficiency 
is  affected  by  changing  the  diffusion  coefficient  and  the  electroos- 
motic  drag  coefficient  of  the  membrane.  Looking  back  to  Eq.  (5), 
we  see  that  these  two  material-based  properties  determine  the 
methanol  concentration  profile  in  the  membrane.  First,  we  set  the 
backing  layer  diffusion  coefficient  to  2.242  x  10-5cm2  s-1  as  in  Eq. 
(20)  and  the  width  of  the  backing  layer  to  0.0180  cm.  Taking  the 
width  of  the  membrane  to  be  0.0179  cm  as  before  and  varying  the 
diffusion  coefficient  from  5  x  10-9cm2  s-1  to  5  x  10-3cm2  s-1,  we 
can  plot  the  fuel  efficiency  as  a  function  of  the  diffusion  coeffi¬ 
cient  and  the  input  feed  concentration  as  in  Fig.  6(a)  and  (b).  As 
expected,  decreasing  the  membrane  diffusion  coefficient,  which  in 
effect  lowers  the  crossover,  has  a  dramatic  positive  effect  on  the 
fuel  efficiency,  especially  in  the  low  current  case.  The  effect  is  more 
pronounced  when  the  load  is  0.05  A  cm-2  because  methanol  is  not 
being  taken  away  as  rapidly  at  the  anode  as  it  is  in  the  0.5  A  cm-2 
case. 

Fig.  6(a)  and  (b)  highlights  again  the  cooperative  effects  of  the 
membrane  properties,  backing  layer  properties  and  system  con¬ 
figuration  on  the  fuel  efficiency.  Altering  the  membrane  diffusion 
coefficient  has  very  little  effect  on  the  fuel  efficiency  except  for  a 
dramatic  change  in  the  region  from  10-4  cm2  s-1  to  10-6  cm2  s-1. 
Any  changes  in  the  fuel  efficiency  outside  this  region  are  limited  by 
the  chosen  backing  layer  properties  and  the  initial  methanol  con¬ 
centration.  The  reader  should  also  note  that  the  data  for  input  feed 
concentrations  of  0.1  M  and  0.5  M  are  not  shown  in  Fig.  4(b).  With 
the  given  set  of  backing  layer  properties  at  a  0.5  A  cm-2  load,  these 
concentrations  result  in  negative  fuel  efficiencies.  In  other  words, 


Fig.  5.  (a)  Maximum  fuel  efficiency  as  a  function  of  limiting  current  with 
/  =  0.05  A  cm-2,  (b)  Maximum  fuel  efficiency  as  a  function  of  limiting  current  with 
7  =  0.5  A  cm-2. 


there  is  not  sufficient  methanol  in  the  feed  stream  to  sustain  the 
desired  current  at  the  anode  across  the  entire  length  of  the  fuel  cell. 
The  same  goes  for  0.1  M  in  Fig.  4(a),  where  there  is  barely  enough 
methanol  reaching  the  anode  to  sustain  the  desired  current,  result¬ 
ing  in  low  efficiency.  Altering  the  membrane  properties  does  not 
change  this  situation,  so  by  this  constraint  we  can  rule  out  using 
these  input  feed  concentrations  to  optimize  the  fuel  efficiency. 

Fixing  the  diffusion  coefficient  at  the  reference  value  of 
4.9  x  1 0-6cm2  s-1  given  in  Eq.  ( 1 9),  we  now  vary  the  methanol  drag 
coefficient  f^e0H  =  £h2o/ch2o  t0  gauge  its  effect  on  the  fuel  effi¬ 
ciency.  For  demonstration  purposes  we  arbitrarily  choose  to  vary 
the  value  of  ^'Me0H  over  a  wide  range*  fr°m  0-1  t°  105  cm3  mol-1. 
Recalling  Fig.  2  should  give  us  a  hint  at  the  results.  Increasing  the 
drag  coefficient  by  a  factor  of  twenty  noticeably  increases  the  extent 
of  methanol  crossover  in  the  membrane,  which  should  decrease 
fuel  efficiency.  Given  that  we  are  varying  the  drag  coefficient  by 
a  factor  of  106,  we  should  see  substantial  changes  in  the  fuel  effi¬ 
ciency. 

However,  Fig.  7(a)  and  (b)  shows  the  maximum  fuel  efficiency 
changes  very  little  until  §Me0H  >  500  cm3  rnoP1  in  the  low  cur¬ 
rent  case  and  §Me0H  >  50  cm3  rnoP1  in  the  high  current  case. 
Measurements  of  the  water  drag  coefficient  in  a  variety  of  DMFC 
membranes  range  from  §h2o  =  2  to  §h2o  =  5  depending  on  tem¬ 
perature  and  the  equivalent  weight  of  the  membrane  [35].  These 
values  correspond  to  a  range  of  §Me0H  =  36  cm3  moP1  to  §Me0H  = 
80  cm3  mol-1,  which  includes  the  reference  value  of  §Me0H  = 
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Fig.  6.  (a)  Maximum  fuel  efficiency  vs.  membrane  diffusion  coefficient  for  varying 
feed  concentrations  and  /= 0.05  A  cm-2,  (b)  Maximum  fuel  efficiency  vs.  membrane 
diffusion  coefficient  for  varying  feed  concentrations  and  /  =  0.5  Acrrr2. 

46  cm3  mol-1  given  in  Section  2.2.  In  other  words,  the  drag  coef¬ 
ficient  or  the  current  must  be  very  high  in  order  for  electroosmotic 
drag  to  affect  efficiency  as  much  as  diffusion  does.  The  expected 
trend  of  decreasing  fuel  efficiency  for  increasing  drag  does  exist, 
but  the  effect  of  drag  is  minimal  within  the  range  of  physically  rea¬ 
sonable  ¥Me0H  values,  especially  in  the  low  current  case.  At  low 
current  densities,  the  given  membrane  and  backing  layer  proper¬ 
ties  must  be  altered  if  the  drag  is  to  have  a  greater  effect  on  the  fuel 
efficiency. 

Throughout  this  analysis  and  especially  in  this  section,  it  is 
important  to  remember  that  this  model  is  concerned  only  with  fuel 
efficiency  and  not  system  efficiency.  System  efficiency  is  a  more 
complete  measure  of  fuel  cell  performance  because  it  is  a  product 
of  the  reversible  thermodynamic  efficiency  (for  an  ideal  fuel  cell), 
the  voltage  efficiency  (incorporating  kinetic  losses),  and  the  fuel 
efficiency.  Optimizing  the  fuel  efficiency  is  not  always  correlated 
to  optimizing  the  voltage  efficiency,  so  achieving  maximum  fuel 
efficiency  does  not  guarantee  maximum  system  efficiency.  Only  in 
applications  where  portability  takes  precedence  over  power  out¬ 
put  is  fuel  efficiency  of  primary  concern,  and  it  is  those  applications 
that  this  model  is  most  useful  in  system  design.  Keeping  this  point 
in  mind,  the  next  section  will  consider  how  to  estimate  the  required 
size  of  system  components  for  a  fuel  cell  operating  at  its  maximum 
fuel  efficiency. 


4.  Methanol  flux  and  the  size  of  system  components 

Because  portability  is  of  primary  importance  in  many  applica¬ 
tions  where  DMFCs  might  be  used,  we  want  to  use  this  model  to 
predict  how  large  the  methanol  storage  tank,  circulation  pump, 


Fig.  7.  (a)  Maximum  fuel  efficiency  vs.  methanol  drag  coefficient  for  varying  feed 
concentrations  and  7  =  0.05  A  cm-2,  (b)  Maximum  fuel  efficiency  vs.  methanol  drag 
coefficient  for  varying  feed  concentrations  and  7  =  0.5  A  cm-2. 

cell  stack,  and  other  components  need  to  be  to  produce  a  required 
methanol  flow  rate.  The  less  methanol  that  is  used  and  the  more 
slowly  it  is  passed  through  the  anode  channel,  the  smaller  we  can 
make  these  various  system  components.  However,  a  minimum  flow 
rate  is  set  by  the  current  that  the  cell  must  generate.  Ideally,  there 
would  be  no  methanol  crossover,  and  all  of  the  fuel  in  the  system 
would  either  be  consumed  in  the  desired  reaction  or  recirculated 
through  the  system  to  be  consumed  later;  the  size  of  the  fuel  tank 
would  be  solely  dependent  on  the  demand  for  current.  However, 
the  existence  of  crossover  and  the  potential  for  wasted  fuel  means 
that  the  size  of  the  tank  also  depends  on  the  fuel  efficiency. 

First,  lets  consider  the  size  of  the  methanol  storage  tank  in  a  situ¬ 
ation  where  the  fuel  is  not  recycled.  The  minimum  flow  rate  into  the 
anode  channel  needed  to  produce  the  required  feed  concentration 
is 

^MeOH  =  C feed  ^channel  V  =  QpS’  (22) 

where  A  is  the  electrode  area,  which  is  equal  to  w  •  Lchannei.  In  this 
case,  the  input  flow  rate  or  flux  is  directly  proportional  to  the  sto¬ 
ichiometric  ratio  and  therefore  inversely  proportional  to  the  fuel 
efficiency.  The  size  of  the  fuel  tank  and  the  circulation  pump  are 
in  turn  proportional  to  the  flow  rate.  Therefore,  we  can  employ 
the  analysis  used  in  Section  3  to  maximize  the  fuel  efficiency  and 
thereby  minimize  the  size  of  the  fuel  tank  and  circulation  pump. 
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Fig.  8.  Schematic  of  the  entire  fuel  cell  system. 


Now  we  consider  a  scenario,  depicted  in  Fig.  8,  where  the 
methanol/water  solution  is  separated  from  the  carbon  dioxide  at 
the  stack  exit  and  is  recirculated  through  the  system.  The  output 
flux  from  the  anode  channel  is  the  input  flux  minus  the  flux  that 
is  redirected  into  the  cell.  The  value  of  co  at  the  stack  exit  repre¬ 
sents  the  fraction  of  methanol  that  is  consumed  in  the  entire  cell, 
so  multiplying  the  input  flux  by  co  gives  us  the  output  flow  rate 

»ss;-jr(»+>*+£)-  1231 

The  required  flow  rate  of  methanol  from  the  storage  tank  into  the 
fuel  system  is  the  input  flow  rate  minus  the  output  flow  rate: 

(24) 

In  this  case,  the  stoichiometric  ratio  is  not  the  only  parameter  that 
determines  the  methanol  flow  rate  and  the  size  of  the  fuel  tank. 
Fig.  9(a)  and  (b)  plots  the  flux  -  without  including  the  area  of  the 
anode  A  as  is  shown  above  -  as  a  function  of  the  limiting  current 
for  a  system  with  recirculation.  Two  points  must  be  considered 
when  looking  at  these  plots.  First,  the  plots  cut  off  when  the  limiting 
current  reaches  the  desired  current  because  allowing  the  limiting 
current  to  reach  or  to  fall  below  the  desired  current  will  produce  a 
nonphysical  result.  Also  note  that  increasing  the  load  requirement 
by  a  factor  of  ten  increases  the  required  methanol  flux  from  the 
storage  tank  by  a  factor  ranging  from  about  five  to  ten.  Designing  a 
fuel  delivery  system  that  can  handle  such  a  wide  range  of  methanol 
flow  rates  is  an  important  engineering  problem  that  this  model 
forecasts  but  of  course  does  not  answer. 

An  important  question  that  this  analysis  can  answer  is  the 
question  of  whether  a  recirculation  system  is  needed.  Comparing 
Fig.  9(a)  and  (b)  to  Fig.  5(a)  and  (b)  demonstrates  the  differences  in 
required  methanol  flux  between  a  system  with  recirculation  and 
a  system  where  the  flux  depends  solely  on  the  efficiency.  In  both 
the  high  and  low  current  scenarios  with  recirculation,  the  flux  - 
and  consequently  the  tank  size  -  increases  as  the  anode  limiting 
current  increases.  Without  recirculation,  the  trend  is  the  opposite. 
Therefore,  one  can  consider  the  tradeoffs  of  a  recirculation  system 
once  one  knows  the  duration  of  operation  and  scaling  of  pump  sizes 
and  tank  sizes  with  total  capacity.  If  the  optimum  conditions  for 
the  system  dictate  that  very  little  fuel  comes  out  of  the  stack,  then 
the  separators  and  pumps  may  be  a  waste  of  important  space  and 
power  in  a  portable  device. 


-2 

Limiting  current/ A  cm 


Fig.  9.  (a)  Storage  tank  flux  as  a  function  of  limiting  current  with  I  =  0.05  A  cm-2  .(b) 
Storage  tank  flux  as  a  function  of  limiting  current  with  I  =  0.5  A  cm-2 . 


5.  Conclusion 

Using  analytic  calculations  of  the  methanol  concentration,  it  is 
possible  to  illustrate  the  effects  of  material  properties  and  system 
design  on  the  fuel  efficiency  of  the  direct  methanol  fuel  cell.  Fuel 
efficiency  can  be  calculated  by  taking  the  ratio  of  the  amount  of 
methanol  that  enters  the  cell  to  the  amount  that  is  consumed  in 
the  desired  reaction.  By  requiring  the  stoichiometric  ratio  to  be 
maintained  above  a  minimum  value  at  the  stack  exit,  a  maximum 
achievable  fuel  efficiency  can  be  calculated.  Optimizing  the  fuel  effi¬ 
ciency  based  on  the  limiting  current  at  the  anode  requires  matching 
the  backing  layer  properties  to  the  feed  concentration  and  to  the 
load  requirement.  Lowering  the  diffusion  coefficient  in  the  mem¬ 
brane  always  improves  the  fuel  efficiency  by  lowering  methanol 
crossover,  especially  in  the  low  current  scenario.  Altering  the  elec- 
troosmotic  drag  coefficient  of  methanol  in  the  membrane  has  very 
little  effect  on  the  fuel  efficiency  unless  either  the  current  or  the 
drag  coefficient  is  very  high.  The  size  of  system  components  can 
also  be  calculated  based  on  the  required  methanol  flux  into  and 
out  of  the  cell  operating  at  maximum  efficiency.  Generally  speak¬ 
ing,  the  size  of  the  fuel  tank  and  other  components  decreases  as  the 
feed  concentration  and  load  requirements  decrease. 

While  this  model  does  not  allow  for  tradeoffs  of  power  versus 
energy  or  the  likely  tradeoff  between  ohmic  losses  and  crossover  in 
the  membrane,  future  work  may  incorporate  electrochemical  equa- 
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tions  into  this  model  so  that  it  can  make  predictions  about  power 
output  and  other  performance  parameters  that  are  important  in 
real  fuel  cell  systems.  Such  an  addition  would  add  complexity  to 
the  model  but  give  it  more  relevance  to  experimental  data.  Despite 
the  model’s  limitations,  it  does  give  insights  into  designing  a  fuel 
cell  system  for  maximum  fuel  efficiency  across  a  dynamic  current 
range.  It  also  offers  pathways  to  more  complex  and  realistic  mod¬ 
eling  that  is  of  interest  to  scientists  seeking  to  maximize  fuel  cell 
performance. 
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Glossary 

i:  current  through  the  external  circuit  (A cm-2) 

ix:  crossover  current  (A cm-2) 

inm:  limiting  current  in  the  backing  layer  (A cm-2) 

ijm;  backing  layer  limiting  current  at  channel  inlet  (Ccm-2  s-1 ) 

\jr.  ratio  of  i°m  to  the  current  through  the  external  circuit  (dimensionless) 
v:  ratio  of  limiting  current  in  backing  layer  and  membrane  (dimensionless) 
Camde •'  methanol  concentration  at  the  anode  (mol cm-3) 

Cfeed:  methanol  concentration  in  the  feed  stream  (mol cm-3) 
c/eed:  feec*  concentration  at  channel  inlet  (z  =  0)  (mol cm-3) 
co:  ratio  of  Cfeed  at  distance  z  compared  to  at  z  =  0  (dimensionless) 

C/,m:  minimum  methanol  concentration  at  the  anode  (mol  cm-3) 

Ch2o-  concentration  of  water  in  the  membrane  (mol  cm-3) 

9:  ratio  of  qim  to  c°eed  (dimensionless) 

D backing-  effective  diffusion  coefficient  in  the  backing  layer  (cm2  s_1 ) 
e:  porosity  of  the  backing  layer  (dimensionless) 

B>MeOH-H2o'-  diffusion  coefficient  of  methanol  in  water  (cm2  s_1 ) 

Dmem:  diffusion  coefficient  in  the  membrane(cm2  s-1 ) 
hacking'-  width  of  the  backing  layer  (cm) 

Lmem:  width  of  the  membrane  (cm) 

NMeOH§:  ^ux  methanol  through  the  backing  layer  (mol  cm-2  s_1 ) 

NMeOH:  ^ux  methanol  through  the  membrane  (mol  cm-2  s-1 ) 

NMeoHel :  ^ux  methanol  through  the  flow  channel  (mol  cm-2  s-1 ) 

£ MeOH •  electroosmotic  drag  factor  of  methanol  (dimensionless) 

$'Me0H:  electroosmotic  drag  coefficient  of  methanol  (cm3  mol-1 ) 

£h2o:  electroosmotic  drag  coefficient  of  water  (cm3  mol-1 ) 

0:  ratio  of  flux  due  to  drag  relative  to  diffusion  (dimensionless) 
n:  mathematical  construct  (no  physical  meaning)  (dimensionless) 
w:  flow  channel  separation  width  (cm) 
z:  distance  down  the  flow  channel  (cm) 

z':  normalized  distance  down  the  flow  channel  (dimensionless) 
s:  stoichiometric  ratio  for  the  methanol  flow  (dimensionless) 

A, channel •  cross-sectional  area  of  the  flow  channel  (cm2) 

A:  superficial  area  of  the  anode  (cm2) 

NMeOH:  methanol  flux  into  the  cell  (mol  cm-2  s_1 ) 

NMeOH:  methanol  flux  out  of  the  cell  (mol  cm"2  s_1 ) 

NMeOH:  methanol  flux  out  of  the  tank  (mol  cm-2  s-1 ) 


